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^ ! Abstract 

' A graviton of a nonzero mass and decay width propagates five physical 

. polarizations. The question of interactions of these polarizations is crucial for 

! viability of models of massive/metastable gravity. This question is addressed 

' in the context of the DGP model of a metastable graviton. First, I argue that 

■ the well-known breakdown of a naive perturbative expansion at a low scale 
'nJ" ! is an artifact of the weak- field expansion itself. Then, I propose a different 

I expansion - the constrained perturbation theory - in which the breakdown 

■ does not occur and the theory is perturbatively tractable all the way up to its 
^ , natural ultraviolet cutoff. In this approach the couplings of the extra polar- 
Qj [ izations to matter and their selfcouplings appear to be suppressed and should 

■ be neglected in measurements at sub-horizon scales. The model reproduces 
J> . results of General Relativity at observable distances with high accuracy, while 

^ I predicting deviations from them at the present-day horizon scale. 



1 Introduction 

An idea that the observable acceleration of the Universe could be a result of large- 
distance modification of gravity is attractive, and is experimentally testable PP- 
0. Moreover, large-distance modifications of gravity give rise to a conceptually 
new approach to the long-standing cosmological constant problem jHl Ej. Hence, 
development of the models of modified gravity becomes an important task. 

The DGP model jHj is a covariant theory of the large-distance modification of 
gravity (see, e.g., [Ol^llSl)- Interactions in this model are mediated by a single 
graviton that lives in infinite-volume five-dimensional space-time. This graviton 
resembles a massive 4D spin-2 state since it has five polarizations. Experimental 
constraints on extra light states with gravitational interactions are rather severe, 
therefore, the question of how these extra polarizations interact with observable 
matter and with themselves becomes crucial for the consistency of the model. 

From a 4D perspective, the graviton behaves as a state of a nonzero decay width^. 
As such it shares certain properties of a massive graviton (and massive non-Abelian 
gauge fields). If non-Abelian gauge fields are prescribed mass my by hand (i.e., 
without using the Higgs mechanism), then the nonlinear amplitudes of the theory 
become strong at a scale my / g-, where (7 is a gauge-coupling constant. A similar effect 
exists in a theory of gravity in which mass rrtg is introduced by hand as in Ref . [2H] ■ 
This was shown in Ref. by calculating nonlinear corrections to a spherically 

symmetric static body of the gravitational radius r^, in which case the weak-field 
approximation breaks down at a scale A^ ~ ^g/ {f'g^gY^^ |2H]- This breakdown is 
a universal property of the perturbative expansion in the massive theory as it can 
be understood by looking at tree- level trilinear graviton vertex diagrams ISO]- A 
simplest trilinear (graviton)^ diagram gives rise to the singularity For 

small irig this diagram is enhanced even though it is multiplied by an extra power 
of the Newton constant. This can also be understood in terms of interactions of the 
longitudinal polarizations becoming strong jST]. For the pure gravitational sector 
itself the corresponding scale A^ reduces to mg/{mg/MpiY^^ , which can be made 
only as big as mg/{mg/MpiY^^ by adding higher nonlinear terms PT]- 

The reason for the breakdown of perturbation theory at a low scale can be traced 
back to terms in the graviton propagator that contain products of the structure 

rrig 

with similar structures or with the flat space metric. These terms do not manifest 
themselves in physical amplitudes at the linear level since they are multiplied by 
conserved currents, however, they enter nonlinear diagrams leading to the breakdown 
of perturbation theory for massive non-Abelian gauge fields or massive gravity ^. 

^In 6D or higher |23, it behaves as a massive graviton with a decay width [251 126[I77| . 

^Unfortunately, the massive gravity in 4D is an unstable theory with an instability 
time scale that can be rather short |33| . 

1 



(1) 



In the model of Ref. [S] the graviton decay width is introduced by adding to 
General Relativity a term that is reparametrization invariant (see the action (^). 
As a result, the terms similar to (jT)) in the propagator are gauge dependent. In a 
simple gauge adopted in Ref. [S] they read as follows: 



where rric ~ 10^'^^ eV is a counterpart of rrig and p = is a root of the Euclidean 
momentum square. Because of the singularity in rric in (j2]), perturbation theory 
breaks down precociously 30j. However, this breakdown is an artifact of an ill- 
defined perturbative expansion - the known exact solutions of the model have no 
trace of the breakdown scale (see Ref. [SO]). This shows that if one sums up all 
the tree-level perturbative diagrams, then the breakdown scale should disappear. 
A more ambitious program is to asks for quantum consistency of the model as an 
effective theory with a cutoff. Such a consistency can only be established within 
the perturbation theory. The latter is ill-defined above the scale mg/ {nig/ M-p\)^/^ 
jHl] (see also [33]), because of the same nonlinear diagrams that make the weak- 
field expansion ill-defined already at the classical level |30j. However, it is not 
clear whether the new scale that emerges in perturbation theory has any physical 
significance. The perturbative breakdown at the quantum level could be an artifact 
of the technical method itself. As was recently argued by Dvali jSH] 5 in a certain toy 
sigma model a similar breakdown takes place, however, the full resumed solution of 
the theory is valid well above the naive breaking scale. 

In the next section I will discuss in detail the reasons for the breakdown of 
perturbation theory in the DGP model. Then, in section 3 I will show how one can 
define a different perturbative procedure such that the singular terms like are 
eliminated from the propagator and the UV behavior of the amplitudes is regular 
all the way up to the cutoff of the theory. The effective one-graviton exchange 
amplitude between two sources with the stress-tensors T^^ that I will calculate in 
section 3 takes the form: 



This amplitude interpolates between the four-dimensional behavior at p ^ rric and 
the five-dimensional behavior at p <^ nic. It has no poles on a physical Riemann 
sheet and satisfies requirements of unitarity, analyticity and causality. 

Notice that the naive perturbation theory in the DGP model has no strong 
coupling scale if the localized 5D Einstein-Hilbert term is used on the brane world- 
volume 1^ and in the "dielectric" regularization of the DGP model jSH] • In higher 
dimensional generalizations of the DGP model [2111221, the perturbation theory is 
well-defined and there is no problem in the first place as was recently shown in 
[SZj. The concern of the present work is the five-dimensional model with the four- 
dimensional induced term (see the action (jH) below). 



rricP 



(2) 





(3) 
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2 What's wrong with perturbative expansion? 



A brief answer to the above question is as follows. The DGP model has two dimen- 
sionful parameters: the Newton constant and the graviton lifetime m^. As a 
result, the naive perturbative expansion in powers of Gtv is contaminated by powers 
of l/rric. Hence, for small values of perturbation theory breaks down for the 
unusually low value of the energy scale. The question is whether this breakdown is 
an artifact of perturbation theory, or it could be that the breaking scale is truly a 
physical scale at which the model needs certain UV completion. Depending on a 
concrete model at hand, the either of above two possibilities could be realized. 
To study these issues in detail we consider the action of the DGP model jS] 

S = Ml, J d'x^R{g) +Ml J d'xdy^n.ig) , (4) 

where R and TZ^ are the four-dimensional and five-dimensional Ricci scalars, respec- 
tively, and stands for the gravitational scale of the bulk theory. The analog of 
the graviton mass is rric = 2M^/Mp,. The higher-dimensional and four- dimensional 
metric tensors are related as 

g{x,y = 0) = g{x) . (5) 

There is a boundary (a brane) at y = and Z2 symmetry across the boundary 
is imposed. The presence of the boundary Gibbons-Hawking term is implied to 
warrant the correct Einstein equations in the bulk. Matter fields are assumed to be 
localized on a brane and at low energies, that we observe, they do not escape into 
the bulk. Hence, the matter action is completely four-dimensional Sm = J d'^xLM- 

Our conventions are: 77^5 = diag [H ] ; A, i? = 0, 1, 2, 3, 5 ; /x, z/ = 0, 1, 2, 3. 

The simplest problem is to calculate the Green's function D^^''^^ and the ampli- 
tude of interaction of two sources T^^ on the brane 

A — 1^ DfJ-'^',°'P (p,\ 

"^1— graviton — fiu -^-^ a/3 • \'J J 

In order to perform perturbative calculations one has to fix a gauge. A simple way 
adopted in JH| is to use harmonic gauge in the bulk d^hAB = dBh(^/2. In this case 
the one-graviton exchange amplitude on the brane (in the momentum space) takes 
the form: 

^1 -graviton (P, y = O) = , (7) 

where we denote the Euclidean four-momentum square by p^: 

2 u 2|2|2|2 2, 2,2,2 /o\ 

p = -p^p^ = -Po + Pi +P2 +P3 = + Pi +P2 +P3^ (8) 

and 

Tl/s = SnG^{T'^^ - It-t) . (9) 
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We can see that ^i-graviton in © is non-singular in the mc — hmit (moreover, the 
pole in the amplitude is on a nonphysical Riemann sheet, see discussions in the next 
section.). However, the gauge dependent part of the propagator Z}^"^'"/^ contains 
terms proportional to ((21) . These terms do not enter the one-graviton exchange 
amplitude, but, they do contribute to higher-order tree- level non-linear diagrams 
which blow up in the limit |3Uj . 

The same fact is reflected in the expression for the trace of h^^ = g^^, — rj^^, which 
in the harmonic gauge takes the form jS]: 

K{p,y = o) = -T7^- (10) 

(Hereafter the tilde-sign denotes Fourier-transformed quantities and we put 8hGn = 
1.). From this expression we learn that: (i) h'^ is a propagating field in this gauge; 

(ii) h'^ propagates as a 5D field, i.e., it does not see the brane kinetic term; (iii) The 
expression for h^^ is singular in the limit rric —>■ 0. The gauge dependent part of the 
momentum-space propagator D{p,y) contains the terms p^puh, which, due to (fTn|l . 
give rise to the singular term Q. Hence, to understand the origin of the breakdown 
of perturbation theory, one should look at the origin of the l/m^ scaling in ()10j). 

The singular behavior of /i^ is a direct consequence of the fact that the four- 
dimensional Ricci curvature R{g) in the linearized approximation is forced to be 
zero by the {55} and/or {fi5} equations of motion. This can be seen by direct 
calculation of R and of those equations, but it is more instructive to see this by 
using the ADM decomposition. The {55} equation reads: 

R=iK:)'-Kl,, (11) 

where K^^, denotes the extrinsic curvature. Since K ~ 0{h) the above equation 
implies that the four-dimensional curvature R ~ 0{h'^) and in the linearized order 
R vanishes. Let us now see how this leads to the singular behavior of h in (fTUI) . The 
junction condition across the brane contains two types of terms: there are terms 
proportional to nic and there are terms that are independent of rric. The former come 
from the bulk Einstein-Hilbert action while the latter appear due to the worldvolume 
Einstein-Hilbert term. In the trace of the junction condition the rric-independent 
term is proportional to the four-dimensional Ricci scalar R. On the other hand, 
as we argued above, R has no linear in h term in the weak-field expansion, simply 
because these terms cancel out due to the {55} and/or {/u5} equations. Therefore, 
in the linearized approximation the junction condition contains only the terms that 
come from the bulk. These terms are proportional to mji. This inevitably leads to 
the trace of h (fTUI) that is singular in the limit and triggers the breakdown 

of the perturbative approach as discussed above. 

The above arguments suggest that the two limiting procedures, first truncating 
the small h expansion and only then taking the tjIc — > limit, do not commute 
with each other. Therefore, the right way to perform the calculations is either to 
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look at exact solutions of classical equations of motion, as was argued in P^ I3Uj. or 
to retain at least quadratic terms in the equations. The obtained results won't be 
singular in the rric limit. 

However, neither of the above approaches addresses the issue of quantum grav- 
itational loops. Since the loops can only be calculated within a well-defined per- 
turbation theory, one needs to construct a new perturbative expansion that would 
make diagrams tractable at short distances. 

In the next section we will propose to rearrange perturbation theory in such a 
way that the consistent answers be obtained in the weak-field approximation. This 
can be achieved if the linearized gauge-fixing terms can play the role similar to the 
nonlinear terms. We will see that this requires a certain nontrivial procedure of 
gauge-fixing and choosing of appropriate boundary conditions. 

3 Constrained perturbative expansion 

Below we develop a perturbative approach that allows to perform calculations in 
the weak-field approximation without breaking the expansion at a low scale. 

We recall that in the DGP model the boundary (the brane) breaks explicitly 
translational invariance in the y direction, as well as the rotational symmetry that 
involves the y coordinate. However, this fact is not reflected in the linearized approxi- 
mation - the linearized theory that follow from Q is invariant under five- dimensional 
reparametrizations^. This line of arguments suggests to introduce constraints in the 
linearized theory that would account for the broken symmetries. It is clear that an 
arbitrary set of such constraint cannot be consistent with equations of motion with 
boundary conditions on the brane and at ?/ — oo. However, by trial and error a 
consistent set of constraints and gauge conditions can be found. Below we introduce 
this set of equations step by step. We start by imposing the following condition: 

= d^h, + d'^K^ = 0. (12) 

Furthermore, to make the kinetic term for the {/u5} component invertible we set a 
second condition: 

B^ = d^h^^ = Q. (13) 

At a first sight, the two conditions (fT^ and (fT^ fix all the x-dependent gauge trans- 
formations and make the gauge kinetic terms non-singular and invertible. However, 
at a closer inspection this does not appear to be satisfactory. One can look at the 
{yuz/} component of the equations of motion and integrate this equation w.r.t. y 
from — e to e, with e ^ 0. After the integration, all the terms with and B^ 

■^If instead of the boundary we consider a dynamical brane of a nonzero tension, then the five- 
dimensional Poincare symmetry is nonlinearly realized and one has to include a Nambu-Goldstone 
mode on the brane. 
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vanish. The resulting equation (which is just the Israel junction condition) taken 
by its own, is invariant under the following four-dimensional transformations 

K'^W^y'>\y=o = V(2^>l/)l3/=o + df,C^\y=o + duC^,\y=o. (14) 

This suggests that in the rric — > limit the gauge kinetic term on the brane is 
not invertible. As a result, the problem of a precocious breakdown of perturbation 
theory discussed in the previous section arises. To avoid this difficulty one can 
introduce the following term on the brane worldvolume: 

AS = -Ml, j d'^xdydiy) [d^h^, - • (15) 

This makes the graviton kinetic term of the brane invertible even in the ttIc — > 
limit. At this stage, the partition function can be defined as follows: 

Zgf = lim„,^^o / dhAB^^Vi{s + AS +Ml j d^xdy + . (16) 

Here S and AS* are given in (0)) and (fT3j) respectively, and the limit a, 7 — >■ enforces 
(fT^ and (fT^. Before proceeding further, notice that Eqs. (fT^ and (fTSj) would have 
been just gauge-fixing conditions if the boundary were absent (e.g., in a pure 5D 
theory with no brane). However, in the present case, the above equations, when 
combined with the junction condition across the brane, enforce certain boundary 
conditions on the brane. Therefore, Eqs. (fT^ and (fT^ do more than gauge-fixing, 
and 7 and a cannot be regarded as gauge fixing parameters. The prescription given 
by (fT^ is to calculate first all Green's functions and then take the limit 0,7 0. 
Because of this, the results of the present calculations differ from (8j where other 
boundary conditions were implied. 

Using (jl6|) we calculate below the propagator D and the amplitude A defined in 
dni)- We will see that there are no terms in D that blow up as nic 0. 

We start with the equations of motion that follow from ()16|). The {/uz/} equation 
on the brane reads 

j dy {dl)h^^ - v^^d^h^ + d^d^h^y + d^d^h^^ - 2 r/^^^^^s/is^ ) + Cjfj 

- {d^dc^hau + dyd^K^ - d^dyh" - ri^yd^dph'^^ + ^Ti^ydjh'^) = T^^, (17) 

where e — 0, (9|, = OaO^, df = dfj,d^. In (|T7j) we retained only terms that are 
nonzero in the e — limit. Furthermore, G^^^ denotes the 4D Einstein tensor: 

eg) = dl V - d^^d^K - d,d^K + K - Vf^^dj K + dp (18) 
The {^iv} equation in the bulk takes the form: 

+ d^d^h^^ - 2 77^^(9" Ss/isc - ^(9^<9^/i55 + df.d'^K^,) = 0. (19) 
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As a next step we turn to the {/i5} equation which reads as follows: 

dlh,^ - d^d^h^, - d,{d^h^, - d^h'-J - - {d.d'^K^) = . (20) 

7 

Finally, the {55} equation takes the form 

dl hi - d^d^h'^'' - -{dl + d,d,hn = . (21) 
a 

After the calculation is done the limit a, 7 — > should be taken. We turn to the 
momentum space w.r.t. four worldvolume coordinates: 

hAB{x,y) = Jd^pe'P^hAB{p,y). (22) 

From the above equations we calculate the response of gravity to the source T^,^. In 
the limit a, 7 — > the results are as follow: 

y) - ^ ~\v,.T 4^1^) e-'^l . (23) 

+ m^p \ 2 p^ + SnicPj 

We note that in this expression there are no terms similar to unlike to what 
happens in the harmonic gauge |H] where the singular terms are present. 
For the off-diagonal components we find that ha^ ~ 7Pa, and 

ha5{P,y)^0- (24) 

Finally, 

rri 

2 " 2p^ + 3m^p ' ^ ^ 

with r = [p"^ + 2mcp) / [p'^ + mcp) . The amplitude on the brane, as was already stated 
in 0, takes the form 



p2 + iTLcP \ ^'^ 2 p"^ + SrricP 



A-graviton(p, V = 0) = ^ - -T' " , , . (26) 



A remarkable property of this amplitude is that it interpolates between the 4D 
behavior at p mc 

AMp,y = 0) - ^ - It') , (27) 



and the 5D amplitude at p <^ rric 
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Figure 1 : The complex plane of Minkowskian momentum square with the branch cut and 
poles on a second Riemann sheet. The physical Riemann sheet is pole free. 

This amplitude has no vDVZ discontinuity |^ |32] ■ 

It is instructive to rewrite the amplitude ()26|) in the following form: 



■^1— graviton 9 , ,^ 9 , ' ) 

+ nicP D + rricP 



where 



and 



= - 1t-T) , (30) 



9{p') - , ""rf ■ (31) 

p'^ + SrricP 

The first term on the r.h.s. of ()29|1 is due to two transverse polarizations of the 
graviton, while the second term is due to an extra scalar polarization. The scalar 
acquires a momentum-dependent form-factor. The form-factor is such that at sub- 
horizon distances, i.e., when p ^ rric, the scalar decouples. At these scales the effects 
of the extra polarization is suppressed by a factor rric/p (e.g., in the Solar system 
this is less than 10~^^). However, the scalar polarization kicks in at super-horizon 
scales, p -C m^, where the five dimensional laws or gravity are restored. 

Let us discuss the above results in more detail. For this we study the pole 
structure of the amplitude (f^ . There are two nontrivial poles 

p^ = —rricP, and = —SrUcp. (32) 



Let us find the positions of these poles on a complex plane of the Minkowskian 



momentum square p^, where = p^exp(— ztt). For this we note that there is a 
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branch cut from zero to plus infinity on the complex plane (see Fig. 1). The pole 
at = is just the origin of the branch cut. Because of the cut the complex 
plane has many sheets (the propagator is multivalued function due to the square 
root in it). It is straightforward to show that both of the poles in (jH^ are on a 
non-physical, second Riemann sheet. Moreover, the positions of these poles are far 
away from the branch cut (usual particle physics resonances appear on non-physical 
sheets close to the branch cut, the above poles, however, are located on a negative 
semi-axis of the second Riemann sheet). Hence, the physical Riemann sheet is pole 
free''. The poles on a nonphysical sheet correspond to metastable states that do 
not appear as in and out states in the S-matrix [l^j- Using the contour of Fig. 
1 that encloses the plane with no poles, and taking into account the jump across 
the cut, the four- dimensional Kallen-Lehmann representation can be written for 
the amplitude ()29|1 . The latter warrants four-dimensional analyticity, causality and 
unitarity of the amplitude (j2SI)^- Although the above interpretation is the only 
correct one, one could certainly adopt the following provisional picture that might 
be convenient for intuitive thinking. The second pole in ()32|) can be interpreted as 
a "metastable ghost" with a momentum-dependent decay width that accompanies 
the fifth polarization and cancels its contributions at short distances. Remarkably, 
this state does not give rise to the usual instabilities because it can only appear in 
intermediate states in Feynman diagrams, but does not appear in the in and out 
states in the S-matrix elements. In this respect, it is more appropriate to think that 
the scalar graviton polarization acquires the form-factor g{p) ()3H1 . 

The above results seem somewhat puzzling from the point of view of the Kaluza- 
Klein (KK) decomposition. Conventional intuition would suggest that the spectrum 
of the KK modes consists of massive spin-2 states. The Kallen-Lehmann represen- 
tation for the amplitude as a sum w.r.t. these massive states would give rise to the 
tensorial structure where the first term on the r.h.s. of (j^^ is proportional to Ty^, 
instead of ^^^^ case, the remaining part of the amplitude on the r.h.s. would 

have a negative sign. This might be thought of as a problem. However, this is not 
so. The crucial difference of the present approach from the conventional KK theories 
is that the effective 4D states are mixed states of an infinite number of tensor and 
scalar modes. What is responsible for the mixing between the different spin states 
is the brane-induced term and the present procedure of imposing the constraints. 

^Because of the i/p^ dependence of the propagator there are two choices of the sign of the square 
root. We choose the sign as above. For this choice the poles do not appear on the physical sheet 
and Euclidean Green's functions decay for large y. However, an opposite choice of the sign of 
can also be adopted. This would correspond to a different branch of the theory. On that branch, 
if we insist on flat brane, we find tachyonic poles on a physical Riemann sheet. This indicates 
that Minkowski space on that branch is unstable. The unstable classical solutions found in Ref. 

do precisely correspond to this choice of the sign of the square root. On that branch one can 
also obtain the selfaccelerated solution without introducing the cosmological constant ^1 ^ . This 
branch is decoupled (at least classically) from the branch that we are discussing in this work. 

priori it is not clear why the theory that is truly higher-dimensional at all scales should have 
respected 4D analyticity and causality. 
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In the covariant gauge that we discuss the trace of h propagates and mixes with 
tensor fields. From the KK point of view this would look as an infinite tower of 
states with wrong kinetic terms. However, at least in the linearized approximation, 
the trace is a gauge artifact (similar to the zeroth component of the gauge field in 
covariantly gauge- fixed QED or QCD). Nevertheless, the effect of the trace part is 
that the true physical eigenmodes do not carry a definite four-dimensional spin of 
a local four-dimensional theory (see also |3Hj)- Because of this there is no reason to 
split the amplitude (pUj) into the term that is proportional to Ty^ and the rest. 

The question of interactions of these states in the full nonlinear theory is not 
addressed in the present work. What happens with the diagrams in which the 
"metastable ghosts" propagate in the loops (the unitarity cuts of which should give 
production of these multiple states) remains unknown. However, since the theory 
possesses 4D reparametrization invariance, we expect that these questions will find 
answers similar to those of non-Abelian gauge fields. Further studies are being 
conducted to understand these issues. 

4 Conclusions 

To summarize briefiy, a new, constrained perturbative expansion was proposed. In 
this approach perturbation theory is well-formulated. The resulting amplitude in- 
terpolates between the 4D behavior at observable distances and 5D behavior at 
super-horizon scales. This is due to the scalar polarization of the graviton that ac- 
quires a momentum-dependent form-factor. As a result, the scalar decouples with 
high accuracy from the observables at sub-horizon distances. 

The model can potentially evade the no-go theorem for massive/metastable grav- 
ity jH], that states that for the cancellation of the extra scalar polarization one 
should introduce a ghost that would give rise to instabilities jH) HH HH] . In the 
present model, at least in the linearized approximation, such instabilities do not 
occur. The convenient (although not precise) picture is to think of a "metastable 
ghost" that exists only as an intermediate state in Feynman diagrams which does 
not appear in the final states at least in the linearized theory. Since this state cannot 
be emitted in physical processes, it does not give rise to the usual instability. The 
latter property is similar to the observation made in the "dielectric" regularization 
of the DGP model in [38 . 

The questions that remain open concern the gauge-fixing and interactions in the 
full non-linear theory where the Faddeev-Popov ghosts are expected to play a crucial 
role. These issues will be addressed elsewhere. 

I would like to thank Gia Dvali and Massimo Porrati for valuable comments, and 
Max Libanov for useful critical remarks. 
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